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ABSTRACT 

We present the first results of the Hubble Space Telescope/Faint Object 
Camera long-slit spectroscopy of the inner 8" of the Narrow Line Region of 
NGC4151 at a spatial resolution of 0".029. The emission gas is characterized 
by an underlying general orderly behaviour, consistent with galactic rotation, 
over which are superposed kinematically distinct and strongly localized 
emission structures. High velocity components shifted up to ~ 1500 kms" 1 
from the systemic velocity are seen, associated with individual clouds located 
preferentially along the edges of the radio knots. Off-nuclear blue continuum 
emission is also observed, associated with the brightest emission line clouds. 
Emission line ratios like [Ne III] A3869/[0 II] A3727, and [O II] A3727/H/3 are 
observed to vary substantially between individual clouds. We advance the 
general picture that, as in other Seyfert galaxies observed with HST (e.g. 
NGC1068, Mrk573), the interaction of the radio jet with the ambient gas 
strongly influences both the morphology and the physical conditions of the 
NLR. 
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1. Introduction 

Ground-based long-slit spectroscopy of the Extended Narrow Line Region (ENLR - a 
few kpc scale) of the Seyfert 1 galaxy NGC4151 has shown that the kinematics, ionization 
structure and morphology of the emission-line gas are consistent with photoionization of 
gas in the galactic disk by an anisotropic radiation field produced by the active nucleus 
(Penston et al. 1990). The situation on the scale of the Narrow Line Region (sub-kpc) 
is however, much less clear as there is a large misalignement between the ENLR at PA ~ 
50° (Perez et al. 1989; Robinson et al. 1994) and the NLR which is co-spatial with the 
radio jet at PA ~ 77° (Pedlar et al. 1993). If both the NLR and ENLR structures are to 
be explained in terms of an "ionization cone" then this misalignment implies either a very 
broad (opening angle ~ 120°) cone of ionizing radiation, intercepting the disk at grazing 
incidence (Pedlar et al. 1992, Vila-Vilaro 1993, Robinson et al. 1994) or a geometry 
which is distinct from the "standard" unified model which has the cone axis fixed to that of 
the radio ejecta (Evans et al. 1993). 

The complex structure of clouds and filaments in the inner 5" of the NLR revealed 
by Hubble Space Telescope (HST) images suggests that interaction with the radio jet may 
have an important influence on the observed morphology. In this Letter, we present the 
first results of new HST long-slit spectroscopy of the central 8" of NGC4151 at a spatial 
resolution of 0".029 which enables us to separate for the first time the emission from the 
individual NLR clouds, and demonstrates that the interaction with the radio ejecta strongly 
affects the local kinematic and ionization conditions of the emission gas. Our results imply 
that it is this interaction, rather than anisotropic illumination, which determines the basic 
structure of the NLR. 

For a distance to NGC4151 of 13.3 Mpc, 0".l corresponds to a linear scale of 6.4 pc in 
the plane of the sky. 

2. Observations and Results 

The spectroscopic data consist of four spectra obtained with the Faint Object Camera 
(FOC) f/48 long-slit facility on board the HST. The slit width is 0".063, and its useful length 
about 10". We have also used FOC f/96 images with filters centered at the [O III] A5007 
emission-line and nearby continuum. Table p] summarizes the relevant information on the 
above data set. The f/96 images were reduced using the standard pipeline procedure, as 
described in the FOC Handbook (Nota et al. 1996). 

The spectra were obtained with the slit at PA = 47°, and cover the interval 3650 - 
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5470 A. The data reduction followed the same procedure as described in Macchetto et al. 
(1997); wavelength resolution is 5.1 A ~ 300 kms" 1 at [O III] A5007, and the spatial PSF 
is of the order of 0".08; flux calibration was obtained using the IUE standard LDS749B. 
By comparing the brightness profile of the [O III] A5007 emission line on the spectra with 
that on the f/96 image, we were able to pinpoint the location of the slit to within 0".06 (2 
pixels at the f/48 resolution). Figure [l] (Plate |l|) shows the slit positions superimposed on 
the [O III] A5007 image. The spectra are located on the nucleus, and 0".23, 0".40 and 0".46 
NW of it, and will be hereafter referred to as POS 1, POS2, POS3 and POS4, respectively. 
The nuclear spectrum was heavily saturated in the central 15 pixels, due to the strong 
continuum emission. The presence of a complex radial velocity structure is evident in 
all the brightest narrow emission lines, such as [O II], [Ne III], and [O III]. To study the 
detailed kinematics of the gas the [O III] A4959,5007 emission lines, which provide the best 
resolution and S/N ratio, were extracted in spatial windows varying from two to 20 pixels 
along the slit, and both lines fitted simultaneously using two to six Gaussian components, 
constraining the physical parameters of the corresponding components of each line to their 
theoretical values (3:1 intensity ratio, 48 A separation, and same FWHM). 

Figure |2| (Plate Q) shows the combined VLA+MERLIN 5 GHz radio map from Pedlar 
et al. (1993), overlaid on the [O III] A5007 image of NGC4151. The registration was made 
assuming that the optical nucleus corresponds to component C4, as advocated by Mundell 
et al. (1995). There is an apparent anti-correlation between the optical and the radio 
emission, with the brightest emission-line filaments surrounding the radio knots as would 
be expected if the radio-jet is clearing a channel in the surrounding medium and enhancing 
the line emission along its edges by compression of the ambient gas. 

3. Kinematics 

Figure shows the results of the [O III] A4959,5007 fitting for three of the four 
slit positions (the POS 3 and POS 4 data are very similar to each other and do not 
give independent information). The rotation curves were constructed assuming that the 
narrowest, closest to the systemic velocity, and usually the brightest component (filled 
squares) represent the "main" NLR gas emission. The curves are plotted as a function 
to the distance to the center of the slit, defined as a line passing through the nucleus, 
perpendicular to the slit direction. The solid line on each panel represents the integrated 
[O III] A5007 spatial profile. 

The general behavior of the main component is consistent with normal galactic rotation 
as defined by the ENLR gas emission (Robinson et al. 1994), with strong but localized 
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perturbations superimposed, represented by one, or sometimes two, narrow (FWHM < 800 
kms -1 ) secondary components either blue- or redshifted relative to the core of the line 
(shown as open circles and squares in Figure ^). In the inner 2" of the nuclear spectrum, the 
behaviour of the emission gas shows a strong deviation from the outermost rotation curve. 

It is important to note that the most kinematically disturbed features are found to 
be associated with the emission clouds located along the edges of the strongest radio 
emission components: the red/blue-shifted knots at ~ 1" SW of the center of the slit in 
the POS2 spectrum are located at one side of the C2 radio component, and are well fitted 
by two Gaussian components centered at +1300 and —800 kms -1 from the core velocity, 
and FWHM ~ 330 and 490 kms -1 , respectively; the high velocity (Av ~ —1700 kms -1 ) 
cloud (b) in POS3 lies also alongside this same radio component; the second-brightest 
radio component C3 is in front of the bright emission complex located 0".2-0".5 NW of 
the nucleus, where we observe a clear double-peaked profile in all emission lines, well 
represented by two narrow (FWHM ~ 430km s _1 ) components, one of which is blueshifted 
by ~ 660 kms -1 from the core velocity. 

All these features are characterized by a very small projected size: the two high-velocity 
clouds in POS2 discussed above are 0".17 to 0".23 (less than 15 pc) across, while the 
double-peaked feature in POS 3 (see Figure f|c) is resolved along a 0".3 (20 pc) bin, giving 
further support to the idea that the NLR morphology is strongly shaped in walls and voids 
by the expanding, fast-moving radio components. 

4. Physical Conditions of the Emission Gas 

If the localized kinematics perturbations described above are a consequence of fast 
shocks created by interactions with ejected radio material the physical conditions in the 
clouds may also be influenced by the shocks. Not only can the shocks compression lead 
to local density enhancements (Taylor, Dyson & Axon 1992) but these may also be 
accompanied by secondary photoionizing radiation (Sutherland, Bicknell & Dopita 1993; 
Morse, Raymond Sz Wilson 1996). 

We clearly see that the emission-line ratios vary substantially on scales of a few 
tenths of an arcsec, indicating that the density and ionization state of the emission gas are 
strongly influenced by the local conditions. We also detect several examples of off-nuclear 
blue- continuum emission, consistently associated with the brightest emission-line knots. 

Figure ^ presents a sample of the spectra observed, corresponding to the regions (a) to 
(c) in Figure |3|. The observed line intensities, relative to H/5 = 1.0, are given in Table |2|. 
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As a reference we included the "average" NLR line intensities, obtained by adding up the 
spatially integrated emission of the four slit positions, and the ground-based data for the 
ENLR (r > 6") from Penston et al. (1990). Errors are of the order of 0.05 or smaller for 
the line ratios, and 0.002 or smaller for the [O III] A4363/[0 III] A5007+4959 ratio. 

Spectrum £|(a) corresponds to the brightest [O III] A5007 knot in POS 2, one of the 
clouds associated with the radio component C3, located at a projected distance of 0".18 
SW (12 pc) of the nucleus. Spectrum |](b) corresponds to the region with the high velocity 
cloud located 1" SW (64 pc) of the center of POS 3. This component can be clearly seen 
in the [O III] profile. The FWHM of the "main" [O III] A5007 emission in this spectrum is 
unresolved (310 kins" 1 ), while spectra from neighboring regions (6 pixels, 0".17 ~ 11 pc 
either side) show a A5007 line twice as broad. The H7/H/3 ratios in these two positions are 
consistent with little or no reddening. All line ratios for Spectrum |](b) are for the core 
component only, and do not include the high velocity cloud contribution. 

Both of these spectra present similar ratios [Ne III] A3869/H/3, while [Ne III]/[0 II] A3727 
in (a) is about twice that in (b). The high value of the A3869/A3727 ratio is a local 
characteristic of region (a) that is not observed in nearby clouds, and there is no clear 
evidence of kinematic perturbations in this spectrum. The fact that He IIA4686 in (a) is 
not significantly stronger than in (b) points against the hypothesis of a matter-bounded 
cloud, which will enhance the [Ne III] emission relative to low-ionization lines (Binette, 
Wilson & Storchi-Bergmann 1996). It is also possible that the interaction with the radio-jet 
would compress the gas, and thus lower the ionizaton parameter, favoring the emission of 
low-ionization lines such as [O II] in region (b). Notice, however, that the A3869/A3727 
ratio for this region is similar to the ratio averaged over all the observed area, where one 
would expect the effects of local density enhancements to have been washed out. 

For these two spectra we also measured ratios R[om] = [O III] A4363/[0 III] A5007+4959 
~ 0.012 and 0.016, respectively. Using the 3-level atom of McCall (1984), for a reasonable 
electron density range of 10 3 to 10 5 cm -3 these values correspond to [O III] temperatures 
Tpm] i n the range 12,000 to 16,000 K, slightly above the limit of T e < 11,000 K predicted by 
most photoionization models, but still compatible with the value found for the purportedly 
photoionized, kinematically undisturbed ENLR. In fact, the overall similarity of the NLR 
line ratios to those of the ENLR suggests that variations in electron denstities created by 
different degrees in shock compression in individual clouds are probably the main cause of 
the observed variations of ionization from knot to knot. 

However, this is not the whole story as shown in Spectrum ^(c) which, in addition to 
double-peaked profiles, shows clear evidence of a localized blue continuum. Such localized 
regions of off-nuclear blue continuum emission are seen in several positions in our spectra, 
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and are always associated to, but not exactly coincident with, an emission-line knot. As we 
mentioned above, this is a characteristic spectral signature of the free-free emission of the 
hot gas created by the fast shocks and therefore provides concrete evidence for the existence 
of a local ionizing continuum which might be playing a significant role in the ionization 
structure of the individual clouds. 



5. Conclusions 

The high spatial resolution long-slit spectroscopy of the nuclear region of NGC4151 
presented in this paper clearly show that the popular picture of anisotropic illumination by 
the nuclear source is overly simplistic to explain the complex morphology and kinematics 
observed in the NLR. Comparison of the radio and optical data shows that the line emission 
is enhanced along the edges of the radio knots, as would be expected to result from the 
interaction between the jet and the surrounding medium, similar to what is observed for 
NGC 1068 (Gallimore et al. 1996; Capetti, Macchetto & Lattanzi 1997; Capetti, Axon & 
Macchetto 1997). 

The general behaviour of the [O III] A5007 velocity field is consistent with disk 
rotation, with local deviations caused by turbulence (probably associated with the 
sweeping/expanding motion of the radio cocoon). We detect several distinct and very 
localized high velocity (-1700 to +1300 kms" 1 ), narrow (FWHM < 600 kms -1 ) emission 
components, both blue and redshifted, confined in projected dimensions of up to 0".3. We 
find no evidence for systematic red/blue-shifted deviations as expected by the presence of a 
generalized outflow, as advocated by Schulz (1990). Detailed modeling of the inner \r\ <^ 4" 
kinematics will be presented in a forthcoming paper. 

The data shows that the emission-line ratios can vary substantially within a few tenths 
of arcsec, indicating that the physical conditions of the emission gas are also strongly 
influenced by the local conditions. Furthermore, we detect off-nuclear blue continuum 
emission, consistently associated with the brightest emission-line knots. 

The clear morphological association between the radio and the optical emission, and 
the highly disturbed kinematical components in the NLR of NGC 4151 revealed by our 
data imply that interaction with the radio plasma rather than illumination shapes the 
NLR. Within the context of the present data, the relative balance between the two possible 
explanations of the ionization structure variations we have discussed, density fluctuations 
and/or local ionization, cannot be resolved as we do not have density measurements of 
the individual clouds. Nevertheless, it is evident that shocks associated with the jet-cloud 
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interactions play an important role in determining the physical conditions of the individual 
emission clouds (Taylor, Dyson & Axon 1992; Sutherland at al. 1993). 

In order to unambiguously determine if local ionization effects play a significant role in 
modifying the ionization structure of individual clouds high quality measurements of the 
electron density in the clouds need to be obtained from future HST spectroscopy. 

We thank the referee for carefully reading the manuscript and making a number of 
constructive comments which improved the clarity of the paper. 
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Fig. 1. — [O III] FOC f/96 image showing the position of the spectrograph slit for the spectra 
listed in Table [I]. The white segments represent the position of the spectra shown in Figure 
[|. The actual length of the slit is longer than represented. North is at the top and East is 
to the left. 

Fig. 2. — Contours of the VLA+MERLIN 5 GHz radio-map and knot identification from 
Pedlar et al. (1993) overlayed on the same image as in Figure |l|. The emission gas tends to 
trace the boundaries of the radio-jet and knots. 

Fig. 3. — [O III] rotation curves for three of the observed positions, superimposed on the 
integrated [O III] A5007 spatial profile (solid line), and plotted as a function of the distance 
to the center of the slit. The main kinematic component, corresponding to gas rotating in 
the plane of disk, is represented by filled squares. The open circles and squares correspond to 
secondary narrow and extra components, respectively (see text for details). Velocity errors 
are estimated from the Gaussian fitting procedure. Notice the steep velocity gradient in the 
inner 2" of the nuclear spectrum. Regions (a) to (c) correspond to the spectra shown in 
Figure [|. The scale to the right refers to the integrated flux in the [O III] A5007 profile, and 
the saturated region in the nuclear spectrum is marked as sat.. 

Fig. 4. — Spectra corresponding to regions (a) to (c) in Figure |3|. Top: region (a), the 
brightest line-emission knot in POS 2, showing strong [O III] A4363 and [Ne III] A3869; 
Middle: region (b), containing the high- velocity knot at 1" SW in the POS 3 spectrum. 
The high velocity blueshifted emission from the knot is clearly visible in both A5007 and 
A4959 lines; Bottom: region (c), revealing double-peaked profiles in all the observed lines. 
The dotted line shows the same spectrum in a five times expanded scale to enhance the 
underlying continuum. 
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Table 1. Journal of Observations 



NGC4151 
Observation 


Date 


HST 

Rootname 


Filter 


Exp. time 

(s) 


Spectra 


1996 Jul 03 


X38I0102T 


F305LP 


1247 






X38I0108T 


F305LP 


697 






X38I0109T 


F305LP 


697 






X38I010AT 


F305LP 


697 


Images 


1994 Dec 03 


X2740101T 


F501N 


2160 






X2740102T 


F4ND+F501N 


900 






X2740108T 


F550M 


1340 
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Table 2. Observed Line Ratios: H/9 = 1.0 



Line 


Spectrum (a) 


Spectrum (b) 


Main 


Spectrum 
Blue 


(c) 

Sum 


NLR a 


ENLR b 


[0 II] A3727 


1.54 


3.33 


2.09 


1.21 


1.79 


2.56 


2.44 


[Ne III] A3869 


1.24 


1.25 


0.91 


0.82 


0.88 


1.07 


1.0 


H 7 


0.44 


0.53 


0.39 


0.34 


0.37 


0.39 


0.42 


[0 III] A4363 


0.19 


0.34 


0.37 


0.27 


0.27 


0.25 


0.24 


He II A4686 


0.14 


0.19 


0.26 


0.31 


0.28 


0.27 


0.38 


[0 III] A5007 


11.9 


10.8 C 


15.2 


10.8 


13.7 


12.5 


12.9 


Ratios 
















A3727/A5007 


0.13 


0.31 


0.14 


0.11 


0.13 


0.21 


0.19 


A3869/A3727 


0.81 


0.38 


0.44 


0.68 


0.49 


0.42 


0.41 


R\o in] d 


0.012 


0.016 


0.014 


0.018 


0.015 


0.015 


0.014 



integrated emission from the four slit positions shown in Figure |T|. 

b From Penston et al. (1990), corresponding to the region 6-20" SW of nucleus. 

c Value does not include the high velocity cloud. 

d R [0m ] = [O III] A4363/[0 III] A5007+4959 
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